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ABSTRACT 
The  very  low  frequency  band  has  been  characterized  by 
its  extremely  reliable  propagation  since  the  early  days  of 
radio.   Recent  development  of  oscillators  accurate  to  a  few 
parts  in  lO-H  have  made  it  possible  to  navigate  a  ship  accu- 
rately by  using  existing  VXF  transmissions.   One  method  of 
such  navigation  is  discussed.   In  particular,  navigation  over 
long  periods  of  time  is  explained,  with  corrections  for  di- 
urnal phase  shift  incorporated  in  the  calculation  of  ships 
position.  Equipment  used  to  obtain  experimental  data  is 
described  and  experimental  results  and  sample  calculations 
are  presented. 
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1.   Introduction. 

The  idea  of  navigating  a  ship  using  very  low  frequency 
transmissions  is  not  new.   Emissions  in  the  VLF  band  are 
particularly  appealing  for  this  use  because  of  their  relia- 
bility and  long  range.  With  recent  development  of  highly 
accurate  and  stable  radio  frequency  oscillators,  VLF  naviga- 
tion has  become  both  possible  and  practical.   For  the  method 
of  this  paper  an  oscillator  is  located  at  the  VLF  transmitter 
to  phase  stabilize  the  master  oscillator.  A  second  oscilla- 
tor is  located  on  the  ship.   Both  of  the  oscillators  are  at 
the  same  frequency  and  the  signal  received  from  the  VLF 
transmitter  is  compared  in  a  phase  comparator  with  the  ship- 
board oscillator.   The  output  of  the  phase  comparator  will 
be  constant  so  long  as  the  ship  does  not  move  toward  or  away 
from  the  transmitter  site.  Any  movement  of  the  ship  will  be- 
come apparent  however,  since  the  phase  of  the  received  signal 
will  shift  relative  to  that  of  the  shipboard  oscillator.  At 
the  output  of  a  phase  comparator  calibrated  in  microseconds, 
this  movement  is  indicated  by  a  change  of  6. IB  microseconds 
for  each  nautical  mile  of  movement  toward  or  away  from  the 
transmitter.   The  accuracy  of  this  method  of  determining  the 
movement  of  the  ship  depends  upon  the  accuracy  of  the  two 
oscillators  and  the  stability  of  the  path  of  propagation.   At 
VLF  frequencies  the  path  of  propagation  is  extremely  stable 
and,  in  general,  does  not  depend  upon  distance  from  the 
transmitter. 

If  the  output  of  the  phase  comparator  at  a  known  ships 
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position  were  noted,  a  line  of  present  position  can  be  ob- 
tained perpendicular  to  the  propagation  path  by  finding  the 
microseconds  change.   If  two  or  more  VLF  stations  are  moni- 
tored, a  fix  on  the  position  of  the  ship  can  be  found  at  the 
intersection  of  these  lines  of  position.  Accuracies  of  less 
than  a  mile  are  possible. 

As  might  be  expected  there  are  many  chances  for  error  in 
this  procedure.   In  many  cases  these  errors  can  be  predicted 
and  compensated  for.   In  others  they  can  be  recognized  and 
the  signal  that  is  in  error  discarded  for  navigational  pur- 
poses until  it  is  once  again  free  of  error. 

In  the  past  few  years  considerable  research  has  been 
done  in  the  area  of  VLF  propagation.  Much  of  this  research 
applies  directly  to  the  VLF  navigation  problem;  specifically 
to  the  prediction  of  changes  in  the  propagation  path.   The 
results  of  some  of  this  research,  as  it  applies  to  ship  navi- 
gation will  be  presented  in  this  paper. 

Blackband'7'  and  Stanbrough(23 )  have  used  this  method  of 
navigation  experimentally  for  various  periods  of  time  with 
good  accuracy.   This  paper  relies  heavily  on  their  work  and 
that  of  others. 


2.   The  Waveguide  Mode  Theory  of  VLF  Propagation. 

Propagation  of  radio  waves  at  VLF  frequencies  is  charac- 
terized in  general  by  long  ranges  with  stable  levels  at  the 
receiver.   Considerable  effort  has  been  expended  in  recent 
years  to  determine  the  mechanics  of  this  propagation.   The 
mathematics  involved  in  explaining  this  have  been  found  to  be 
most  complex.   This  complexity  is  primarily  caused  by  the 
curvature  of  the  earth  and  the  non-homogeneous  character  of 
the  reflecting  surfaces  of  the  earth  and  ionosphere.  Many 
papers  have  been  presented  that  take  these  factors  into  con- 
sideration^) (3)  .   In  these  papers  the  volume  between  the 
surface  of  the  earth  and  the  ionosphere  is  considered  to  act 
as  a  waveguide  at  very  low  frequencies.   Some  of  the  apparent 
inconsistencies  of  propagation  at  VLF  are  most  easily  ex- 
plained by  this  waveguide  theory.   Values  of  phase  variations 
and  signal  intensity  at  various  distances  calculated  using 
the  waveguide  model  agree  very  closely  with  experimental  ob- 
servation(lO) (^' •  No  attempt  will  be  made  in  this  paper  to 
justify  the  waveguide  mode  theory  of  VLF  propagation.  Many 
of  the  books  and  papers  listed  in  the  Bibliography  contain 
considerable  mathematical  development  related  to  the  theory. 
The  waveguide  mode  theory  explains  the  characteristics  of  VLF 
propagation  that  are  of  particular  importance  in  the  predic- 
tion of  phase  changes  at  the  receiver.  The  predictions  of 
phase  changes  are  of  great  importance  in  VLF  navigation,  and 
in  this  regard  the  results  of  the  waveguide  mode  theory  will 
be  used. 


3.   The  Diurnal  Phase  Change. 

One  of  the  primary  characteristics  noted  in  VLF  propaga- 
tion studies  is  that  of  the  diurnal  phase  change.   The  propa- 
gation path  between  transmitter  and  receiver  is  extremely 
stable  while  the  entire  path  is  in  daylight.  The  path  for  a 
night  time  propagation  path  is  considerably  less  stable,  and 
in  addition  is  characterized  by  a  phase  shift  from  the  day- 
light path.  The  magnitude  of  this  phase  change  depends  upon 
the  distance  between  the  transmitter  and  the  receiver.  An 
idealized  example  of  a  24  hour  plot  of  the  phase  at  a  re- 
ceiver is  shown  in  Figure  1  (a).  This  shows  an  east  to  west 
transmission  path,  that  is  the  receiver  is  west  of  the  trans- 
mitter. E  to  A  indicates  the  phase  at  the  receiver  during  a 
period  when  the  entire  propagation  path  is  in  daylight. 
Point  A  represents  sunset  at  the  transmitter  site.  The  phase 
changes  linearly  until  B  which  represents  sunset  at  the  re- 
ceiver site.  Line  B-C  represents  the  phase  of  the  received 
signal  during  the  time  the  propagation  path  is  in  darkness. 
Point  C  represents  sunrise  at  the  transmitter  and  the  phase 
changes  to  point  D  which  is  the  time  of  sunrise  at  the  re- 
ceiver. The  phase  of  the  received  signal  is  now  identified 
with  its  previous  value  and  maintains  its  daylight  value  to 
Point  E.  During  the  next  24  hour  period  the  pattern  will  re- 
peat. Figure  1  (b)  shows  one  of  the  many  variations  in  this 
pattern.   In  this  case  the  location  of  the  transmitter  is 
north  and  east  of  the  receiver  so  that  sunset  arrives  at  both 
points  simultaneously. 
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FIG.   1.   DIURNAL  VARIATION  PATTERNS 
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FIG.      1.      DIURNAL   VARIATION  PATTERNS 


The  diurnal  shift  can  be  explained  by  the  waveguide 
theory.   During  the  day,  the  action  of  sunlight  forms  an  ion- 
ized layer  at  an  altitude  of  about  70  kilometers.   During 
daylight  the  walls  of  the  waveguide  are  formed  by  this  layer 
and  the  surface  of  the  earth.   At  night  this  ionized  layer 
disappears  and  the  nighttime  reflecting  surface  appears  at  an 
altitude  of  about  90  kilometers.  When  the  walls  of  the  wave- 
guide are  brought  closer  together,  the  phase  velocity  of  the 
waves  propagating  in  the  guide  is  increased.   Consequently 
the  phase  velocity  of  the  VLF  waves  is  greater  through  a  day- 
light propagation  path.   It  is  convenient  to  refer  to  the 
change  in  phase  as  a  change  in  transmission  time.   This  is 
indicated  as  At   in  Figure  1. 

When  VLF  transmissions  are  used  for  ship  navigation, 
these  diurnal  phase  shifts  must  be  anticipated  and  their  time 
of  occurence  must  be  predicted  so  that  measurements  taken 
during  the  diurnal  shift  may  be  properly  processed.   If  it  is 
desired  to  use  both  day  and  night  transmissions  for  naviga- 
tion, the  pattern  and  magnitude  of  the  shift  must  be  pre- 
dicted in  addition  to  the  times  of  occurrence. 


4.   Prediction  of  the  Diurnal  Change. 

For  purposes  of  prediction  of  the  diurnal  phase  change, 
the  great  circle  path  around  the  earth  may  be  divided  into 
three  general  zones  of  similar  propagation  conditions.   These 
are  shown  in  Figure  2  and  can  be  named: 

Zone  I   -  The  multi-mode  zone 

Zone  II  -  The  single  mode  zone 

Zone  III  -  The  multi-path  zone 
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Figure  2,  Zones  of  similar  propagation  conditions 


4a.  The  Multi-Mode  Zone. 

The  height  of  the  earth-ionosphere  waveguide  is  about  2 
to  9  wavelengths  long  in  the  frequency  range  10  to  30  kilo- 
cycles.  Because  of  this,  many  modes  of  high  order  can  be 
propagated.  These  high  order  modes  are  attenuated  more  rap- 
idly than  the  first-order  mode  however,  and  they  become  small 
relative  to  the  first-order  mode.   The  distance  where  they 
may  be  neglected  depends  upon  the  height  of  the  ionosphere 
reflecting  surface  and  upon  the  frequency  of  transmission, 
but  it  is  generally  considered  to  be  in  the  neighborhood  of 
4000  kilometers'^) •  Within  this  zone  the  modes  interact  in 
a  complex  manner  and  the  phase  of  the  received  signal  may 
vary  in  a  similarly  complex  manner ' 13 )„  The  use  of  VLF  emis- 
sions within  the  multi-mode  region  for  navigation  of  a  moving 
ship  does  not  seem  promising  at  the  present  time  because  of 
this  interaction  between  modes.   Sufficient  VLF  signals  are 
available  from  outside  any  particular  zone  however,  so  this 
is  not  considered  to  be  a  great  disadvantage.   It  is  also 
possible  that  further  research  on  propagation  within  the 
multi-mode  region  will  allow  some  parts  to  be  used.  Caution 
must  be  observed  when  using  VLF  signal  for  navigation.   If 
the  ships  track  passes  within  the  multi-mode  zone  specified 
for  the  VLF  transmitter,  the  signals  received  will  not  be 
usable  for  navigation  in  the  same  manner  as  those  outside 
the  zone. 
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4b.  The  Single  Mode  Zone. 

Beyond  the  multi-mode  zone  the  higher  order  modes  are  no 
longer  of  sufficient  magnitude,  and  only  the  first-order  mode 
need  be  considered.   The  single  mode  region  extends  through- 
out both  zones  II  and  III  in  Figure  2.   This  region  is  char- 
acterized by  VLF  propagation  that  is  normally  stable  with  a 
predictable  diurnal  phase  shift.  This  diurnal  shift  depends 
upon  the  frequency  and  the  distance  from  the  transmitter  to 
the  receiver.  Wait'-^'  gives  an  arithmetical  expression  for 
/D,  derived  from  the  mode  theory; 
^t  =  £p-r  (^  -  £, )  microseconds 
where  D  is  the  distance  to  the  transmitter 

vn  and  vd  are  the  night  and  day 

phase  velocities  respectively. 

c  is  the  speed  of  light  in  space. 
A  value  of  phase  velocity  at  the  frequency  of  interest  can  be 
calculated  using  the  mode  theory,  and  Waitv  '  has  published 
some  charts  of  calculated  phase  velocities  vs.  frequency  for 
various  reflecting  surface  heights.   Figure  3  is  a  plot  of 
these  calculated  values  for  a  perfectly  conducting  ground  and 
imperfectly  conducting  ionosphere.  This  is  a  close  approxi- 
mation to  propagation  over  sea  water.  Measurements  of  the 
phase  velocity  of  VLF  propagation  have  been  experimentally 
measured (22'  at  a  frequency  of  18  kc  and  vd/c  and  vn/c  were 
found  to  be  0.99#  and  0.995  respectively.  These  are  very 
close  to  Wait's  calculated  values  for  h=72  kilometers  during 
daytime  and  h=96  kilometers  during  the  night.   In  the  single 
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FIG.   3.   PHASE  VELOCITY  OF  THE  FIRST  MODE  FOR  A  PERFECTLY 
CONDUCTING  EARTH  AND  AN  IMPERFECTLY  CONDUCTING 
IONOSPHERE  (FROM  WAIT) 
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mode  zone  the  diurnal  phase  shift  increases  nearly  linearly 
with  distance,  and  the  magnitude  of  the  shift  can  be  calcu- 
lated from  At  =  Dkv  where  D  is  the  distance  from  the  trans- 
mitter and  kv  is  a  constant  that  will  depend  upon  the  fre- 
quency.  Table  I  gives  values  of  kv  for  some  frequencies  of 
interest,  where  hd=72  and  hn=96,  kilometers.   These  are  cal- 
culated from  Figure  3. 


Frequency 

£  _  2. 
vn  " '  vd 

kv~ 
(xlO3) 

16.0 

0.00231 

7.70 

17.2 

0.00232 

7.74 

IS. 6 

0.00247 

7.90 

20.0 

0.00241 

S.05 

21.4 

0.00242 

2. 06 

24.0 

0.00262 

S.74 

26.1 

0.00253 

S.60 

TABLE  I 
kv  for  1^=72  and  hn=96  Kilometers 

Blackband  has  determined  experimentally'  1°)  and  V/ait  has 
verified  mathematically'-^)  the  form  of  the  variation  of  At 
vs.  D.  Figure  4  is  a  plot  of  Blackband' s  experimental  data 
at  16  kc.   Blackband  has  chosen  a  linearizing  factor,  kf,  to 
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FIG.     4.%   At   VS,    D  AT  16  KCS. 
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show  the  variation  of  r*t .  Wait*   '  has  found  values  for  kf 
mathematically  from  mode  theory.   These  values  agree  closely 
with  Blackbands'  experimental  kf.   ky  is  related  to  kf  by 
the  relation: 

ky  =  !  7.37  )  x  10-3 


Table  II  gives  values  of  kv  found  from  Wait's  calculated 
values  of  k~. 


f 

(kc) 

kf 

(xl<$3) 

10 

0.7 

11.25 

12 

0.9 

3.75 

14 

1.0 

7.37 

16 

1.03 

7.23 

13 

1.09 

7.22 

13.6 

1.03 

7.23 

20 

1.07 

7.35 

21.4 

1.03 

7.65 

24.0 

0.96 

3.20 

26.1 

0.92 

3.55 

TABLE  II 
Values  of  kv  Calculated  from  kf 
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The  values  of  ky  in  Table  II  seem  to  be  closer  to  those  ob- 
tained experimentally.  Unfortunately  there  are  few  data 
published  of  measured  values  of  A  t  above  1$  kc.   Most  data 
from  experimental  measurements  seem  to  be  for  frequencies  of 
16  and  18  kc.  Even  here,  large  discrepancies  for  ky  appear* 
Tables  III  and  IV  compare  the  values  of  ky  obtained  at  these 
frequencies  at  various  distances.   The  data  is  from  Black- 
band(lg).   Researchers  at  the  ROYAL  AIRCRAFT  ESTABLISHMENT 
(FARNBOROUGH)  have  analyzed  these  and  other  data  to  determine 
if  there  is  any  predictable  effect  on  A  t  caused  by  latitude 
or  direction  of  path.   No  such  effect  has  been  found,  and  At 
appears  to  be  independent  of  these  factors.   Until  more  data 
is  accumulated  on  the  variation  of  At  with  frequency,  it 
appears  that  an  accurata  prediction  of  the  diurnal  phase 
shift  at  a  random  location  is  not  yet  possible • 

At  a  particular  location  however,  the  diurnal  phase 
shift  is  remarkably  constant  from  day  to  day.   Figures  5,  6, 
and  7  are  diurnal  phase  shift  patterns  of  NBA,  NSS  and  NPM 
recorded  at  Monterey,  California  over  several  days.   These 
have  been  corrected  for  frequency  drift.   Figure  8  is  an  ex- 
ample of  an  uncorrected  pattern.   These  are  typical  of  the 
patterns  recorded,  and  the  pattern  of  each  set  is  clear. 
Times  of  sunrise  and  sunset  can  easily  be  seen  on  the  fig- 
ures. These  constant  and  predictable  patterns  of  day  to  day 
diurnal  shifts  have  been  noted  by  a  number  of  researchers. 

Some  of  the  values  of  At  measured  at  Monterey,  Califor- 
nia are  shown  in  Table  V,  and  compared  here  with  the  values 
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that  are  predicted  from  the  value  of  kv  from  Table  II. 


D 
(km) 

At 

(«  sec) 

kv  , 
(xlO^) 

5174 

34o5 

6.66 

5773 

43 

7.45 

6947 

50 

7020 

7424 

43 

6o47 

10901 

72 

6o60 

TABLE  III 


Measured  Values  of  kv  at  16  kcs 


D 
(km) 

At 

(  ju  sec) 

(xlO3) 

4297 

29 

6»91 

4670 

27 

5o3 

3462 

66 

5o7 

3675 

70 

6.07 

9037 

75 

6.30 

9795 

71 

6.95 

TABLE  IV 


Measured  Values  of  kv  at  13  kcs 
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Station 

Frequency 
fkc) 

Distance 
(km) 

Measured 
L   t 

Computed 
At 

GBR 

16 

3591 

52 

62  c  5 

NSS 

21c4 

3970 

40 

30 

NBA 

24 

5242 

42 

43 

NPM 

26ol 

3B63 

43 

33 

TABLE  V 
Comparison  of  Some  Values  of  A  t  Measured  at  Monterey 

With  the  Computed  Values 


The  measured  values  of  A  t  for  both  NSS  and  NPM  are  both 
about  130  percent  of  the  computed  value „   In  these  cases  it 
may  be  that  the  distance  is  not  sufficient  for  the  higher 
order  modes  to  be  attenuated  sufficiently,  and  Monterey  may 
be  within  the  interference  region  for  these  stations,. 

Once  the  diurnal  pattern  for  an  area  is  found,  this 
pattern  can  be  used  to  correct  the  microseconds  readings  for 
that  area.  The  pattern  must  be  modified  as  the  times  of  sun- 
set and  sunrise  change,  and  possibly  slight  modifications  in 
At  will  be  necessary  as  the  seasons  change 0   A  single  pat- 
tern  may  be  used  in  an  area  of  several  hundred  square  miles, 
and  slight  modifications  of  ^t  with  distance  are  possible 
to  extend  a  pattern  to  adjacent  areas . 
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4c  The  Multi-Path  Zone. 

Zone  III  in  Figure  2  is  a  region  where  particular  care 
must  be  taken  when  using  a  VLF  signal  for  navigation.   This 
is  because  of  the  possibility  of  receiving  and  tracking  a 
signal  around  the  longer  great-circle  path.   The  VLF  signal 
around  the  long  great  circle  path  is  always  present  of 
course.   By  using  time  sharing  and  directional  antennas, 
Westfall^l)  in  San  Diego  has  measured  and  compared  phase  and 
amplitude  of  the  signal  transmitted  from  NAA  over  compli- 
mentary great  circle  paths.   Normally  however,  the  signal  re- 
ceived over  the  shorter  path  is  the  greater  in  amplitude  and 
overrides  the  weaker  signal.   In  Zone  III,  the  long  path  may 
be  the  preferred  path  for  various  reasons. 

It  has  been  shown C 17 ) (16)  that  VLF  waves  are  attenuated 
less  when  traveling  from  West  to  East.   Because  of  this  the 
two  path  region  is  tilted  as  shown  in  Figure  1.   In  addi- 
tion, the  boundary  is  not  symmetrical,  and  depends  somewhat 
on  latitude,  the  differences  in  attenuation  being  largest  for 
an  equatorial  path' 3), 

Another  factor  to  be  considered  is  the  difference  in 
phase  velocities  between  nighttime  and  daylight  paths.   Since 
the  phase  velocity  of  the  daylight  path  is  greater,  this  path 
would  be  preferred  over  the  night  path. 

Measurements  have  been  made  of  the  NPM  signal  received 
at  Singapore ( 14)  over  a  period  of  two  years.   The  diurnal 
variation  shows  a  double-humped  pattern  that  is  the  result  of 
splitting  the  propagation  path  into  approximately  two  12= 
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hour  periods.  This  indicates  that,  in  this  case,  the  long 
path  of  29,000  kilometers  is  preferred  over  the  short  path 
of  11,000  kilometers  during  the  daylight. 

The  usefulness  of  VLF  signals  in  Zone  III  for  naviga- 
tional purposes  is  marginal  at  the  present  time,  since  the 
use  of  directional  antennas  on  ships  is  impractical.   Singa- 
pore-NPM  readings  did  exhibit  a  regular  pattern  however,  and 
it  may  be  that  the  propagation  path  may  be  predictable  and 
consequently  a  signal  received  in  Zone  III  usable  for  naviga- 
tion over  long  periods  of  time.   Further  data  appears  neces- 
sary to  determine  this,  and  also  to  determine  the  boundaries 
of  Zone  III  for  each  VLF  transmitter. 
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5o   Seasonal  Variations  of  the  Diurnal  Phase  Change. 

An  obvious  variation  of  the  diurnal  phase  change  with 
the  season  of  the  year  is  that  change  associated  with  the 
length  of  the  day.   During  the  summer  the  daylight  path  is  of 
considerably  longer  duration  than  the  daylight  path  in  the 
winter. 

Observation  of  the  phase  variations  over  various  paths 
indicate  there  are  probably  small  variations  of  the  diurnal 
shift  with  the  change  in  season^),  particularly  for  paths 
that  pass  near  the  polar  regions.   The  height  of  the  iono- 
sphere appears  to  change  with  the  seasons  causing  a  change  in 
phase  velocity.   Researchers  in  Paris,  France  have  measured 
the  diurnal  shift  of  NBA  on  1&  kcs  over  a  period  of  three 
years.  They  have  noted  At  to  be  70  A.  sec  at  the  winter  sol- 
stice, 80  ^Usec   at  the  summer  solstice  and  75  ^sec  at  the 
equinoxes (1° ' •   As  information  on  the  shift  in  phase  velocity 
with  the  change  in  seasons  becomes  available,  a  small  correc- 
tion for  this  obtained  from  published  tables  can  be  included 
in  the  calculations  of  the  line  of  position. 

A  second  and  even  smaller  variation  may  be  caused  by 
changes  through  the  eleven  year  sunspot  cycle.   Sufficient 
information  is  not  available  on  this  point,  but  a  correction 
for  this  can  also  be  easily  made  from  tabulated  values. 
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6.   Anomalies  of  Propagation. 

There  are  several  events  that  may  disrupt  the  path  of 
propagation  at  very  low  frequencies.   These  events  will  con- 
sequently cause  errors,  possibly  of  large  magnitude,  in  the 
VLF  navigation  system.   The  events  are  called  sudden  iono- 
spheric disturbances  (SID).  There  are  three  major  causes  of 
a  SID.   These  are  solar  flares,  magnetic  storms  and  high  al- 
titude nuclear  bursts.   The  effects  of  all  of  these  are  simi- 
lar in  that  they  perturb  the  ionosphere  and  cause  unpredict- 
able discontinuities  of  the  phase  of  the  VLF  signal  received. 
The  length  of  this  perturbation  may  be  from  a  few  hours  for 
a  solar  flare  to  a  few  days  for  a  magnetic  storm.  During 
high  altitude  atomic  tests,  some  paths  have  been  disturbed  as 
much  as  twelve  days,  with  the  height  of  the  ionosphere  de- 
creased by  as  much  as  5  km  from  its  normal  altitude '  *■' )  t      it 
is  not  necessary  that  the  disturbance  occur  directly  along 
the  path  of  propagation  because  of  the  waveguide  nature  of 
the  propagation  involved' 19) 0 

During  these  periods  of  disturbance,  VLF  signals  are  un- 
usable for  navigation,  although  they  are  still  usable  for 
communications.  When  such  a  disturbance  occurs,  ships  using 
such  a  VLF  navigation  system  must  be  immediately  notified  of 
the  disturbance  and  notified  once  again  when  the  ionosphere 
has  returned  to  normal. 

The  frequency  of  occurrence  of  solar-caused  sudden  iono- 
spheric disturbances  of  sufficient  intensity  to  cause  phase 
errors  in  VLF  transmissions  does  not  appear  to  be  known. 
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They  will  occur  most  frequently  during  periods  of  high  solar 
activity.   In  general  it  seems  they  will  occur  infrequently 
and  will  not  be  a  major  source  of  error  to  the  navigation 
system. 
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7.   The  Frequency  Standard,, 

The  frequency  standard  is  the  heart  of  the  VLF  naviga- 
tion system,  and  it  is  only  with  recent  development  of  highly 
accurate  frequency  standards  that  such  a  system  has  become 
possible.   The  standard  used  should  be  of  high  quality  and 
in  any  case  should  have  a  long  term  stability  of  at  least  a 
few  parts  in  10~10.   At  the  U.  S.  Naval  Postgraduate  School 
a  VARIAN  type  V-4700  Rubidium  Frequency  Standard  was  used, 
and  a  standard  of  this  type  seems  well  suited  for  shipboard 
use.   This  particular  standard  is  specified  to  maintain  a 
long  term  stability  of  5x10"  ^   in  a  one  year  period  (stand- 
ard deviation) • 

Whatever  frequency  standard  is  used,  the  initial  pro- 
cedures to  ready  it  for  use  are  similar.   Before  the  stand- 
ard can  be  used  for  navigation  it  must  be  energized  for  a 
sufficient  period  of  time  for  it  to  stabilize.   This  may  be 
a  few  months  for  a  crystal  oscillator  standard  or  only  a  few 
days  for  a  standard  utilizing  an  atomic  resonance  stabiliza- 
tion similar  to  the  Rubidium  Vapor  standard. 

A  short  failure  of  the  power  supplied  to  the  frequency 
standard  will  nullify  the  previous  data  obtained  and  will 
require  that  a  new  ships  position  be  obtained  by  other  means. 
A  failure  of  sufficient  duration  to  allow  the  ovens  to  cool 
will  cause  instability,  and  the  standard  will  not  be  usable 
until  it  has  been  energized  a  sufficient  time  for  it  to 
stabilize  once  again.   Because  of  this,  it  is  standard  prac- 
tice to  supply  the  frequency  standard  from  a  wet  battery 
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supply  and  to  constantly  charge  these  batteries  from  the 
local  source  of  power*   In  the  event  of  a  power  failure,  the 
batteries  supply  full  power  to  the  standard  until  the  main 
supply  is  restored*, 

Once  installed  on  a  ship,  the  frequency  standard  is 
useful  for  other  purposes,,   Radio  transmitters  or  receivers 
can  obtain  their  frequencies  from  a  frequency  synthesizer 
energized  from  the  ships  frequency  standard*   As  an  example 
of  the  accuracy  possible,  a  transmitter  and  receiver  could  be 
set  at  20  megacycles  and  would  be  as  accurate  as 
20,000,000.00  cycles  per  second,  +  0o001  cycle,.   Electric 
clocks  can  be  powered  by  a  source  excited  by  the  frequency 
standard.   Once  the  clocks  are  set  they  would  gain  or  lose 
about  10  microseconds  a  day,  or  a  second  every  300  years 0 
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&.   The  Tracking  Receiver. 

The  second  unit  used  for  obtaining  the  VLF  navigation 
fix  is  the  tracking  receiver.   This  is  a  special  VLF  re- 
ceiver that  measures  the  phase  difference  between  the  local 
oscillator  and  the  signal  received.  At  the  U.  S.  Naval 
Postgraduate  School  a  TRACOR  model  599G  VLF  receiver  was 
used.   This  receiver  was  developed  for  the  purpose  of  uti- 
lizing the  frequency  stabilized  VLF  stations  for  frequency- 
calibration,  and  it  is  particularly  adaptable  for  use  in  the 
VLF  navigation  problem.   A  useful  feature  of  this  receiver 
is  the  digital  dial  on  the  front  panel  that  displays  contin- 
uously the  relative  time  difference  between  the  local  fre- 
quency standard  and  the  incoming  signal,  recording  changes  in 
phase  as  small  as  0.1  microseconds.   Outputs  are  also  avail- 
able to  record  phase  difference  on  an  external  chart.   The 
tracking  receiver  also  electronically  phase  shifts  the  ref- 
erence signal  within  the  receiver  so  it  is  phase  locked  to 
the  incoming  signal.   Figure  9  is  a  simplified  block  diagram 
of  a  tracking  receiver.   Since  the  receiver  bandwidth  is 
narrowed  by  the  phase  tracking  servo,  extremely  high  sensi- 
tivity is  possible.  The  TRACOR  receiver  is  rated  at  a  sensi- 
tivity of  0.01  microvolt  and  will  phase  track  with  an  input 
signal  to  noise  ratio  of  -50db. 

NPM  and  NAA  present  special  problems  during  the  periods 
they  transmit  with  frequency  shift  keying.   It  is  possible 
to  track  a  VLF  station  utilizing  FSK  if  one  side  of  the 
keying  signal  is  locked  to  the  frequency  stabilized  carrier, 
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Experience  at  the  U„  S.  Naval  Postgraduate  School  indicated 
that  NPM  could  be  tracked  during  the  FSK  transmissions  and 
that  NAA  could  not  be  tracked  during  FSK.   Figure  10  is  a 
portion  of  a  chart  of  NPM  received  at  the  U.  S,  Naval  Post- 
graduate School.   NPM  was  transmitting  with  frequency  keying 
for  one  half  hour  commencing  every  odd  hour  (Monterey  time). 
The  FSK  circuitry  caused  an  offset  of  about  10  microseconds 
but  the  signal  was  tracked  and  could  be  used  for  navigation. 

It  will  be  necessary  to  have  a  receiver  for  each  station 
tracked  for  VLF  navigation,,  A  single  frequency  standard  is 
sufficient  for  all  the  receivers.  Normally  the  output  of 
each  receiver  will  be  plotted  on  a  recorder,  and  this  record 
will  indicate  the  microseconds  change  in  phase  as  the  ship 
moves  relative  to  the  VLF  transmitter.  The  microseconds  dial 
on  the  receiver  front  panel  will  also  indicate  this  change. 
The  chart  will  be  useful  to  indicate  the  onset  of  the  diurnal 
phase  change  however. 

A  tracking  receiver  will  give  a  phase  error  that  may  be 
used  for  navigation  from  any  frequency  stabilized  VLF  CW 
transmission.  The  long  time  constant  used  will  completely 
smooth  out  any  CW  keying  the  station  may  be  using.  The 
transmissions  of  most  of  the  VLF  stations  presently  trans- 
mitting can  therefore  be  used  without  change.  Table  VI  lists 
the  VLF  stations  that  are  now  frequency  stabilized.   The  fre- 
quencies of  the  Navy  transmitters  are  subject  to  change. 
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fig,  10.  a  Phase  plot  of  npm  ( 26. ikcs)  while  transmitting 
both  cw  and  fsk  ■ 
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1 
Station 

Freq. 
(kc) 

Location 
(Lat.  &  Long. ) 

Sponsor 

Type  of 
Keying 

WWVL 

20.0 

(40  40!55.9"W) 
(105  02*53.9"W) 

National 
Bureau 
of  Stand. 

CW 

NBA 

2/foO 

(9   04'30"N) 
(79  34f30"N) 

USN 

CW 

NPM 

26.1 

(21   25*N) 
(15S  90'W) 

USN 

CW 
FSK 

NAA 

17.3 

(44  39'N) 
(67  17  »W) 

USN 

CW 
FSK 

NLK 

IS. 6 

(4S   12'N) 
(121  55?W) 

USN 

CW 

NSS 

21.4 

(33  59'N) 
(76  27?W) 

USN 

CW 

GBR 

16.0 

(52  22»N) 
(1   ll'W) 

British 

Post 

Office 

CW 

TABLE  VI 


Frequency-Stabilized  VLF  Stations 
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9.   Determining  Oscillator  Drift  Rates. 

After  the  shipboard  oscillator  has  stabilized,  the  rate 
of  drift  between  this  oscillator  and  the  master  oscillators 
at  the  various  transmitters  must  be  determined.  This  rate  of 
drift  will  vary  from  station  to  station,  but  should  be  steady 
and  predictable.   It  is  this  feature  of  a  small,  predictable 
drift  between  oscillators  that  allows  stabilized  VLF  trans- 
missions to  be  used  for  navigation. 

Pierce  was  the  first  to  discover  that  oscillators  could 
be  compared  with  high  precision  through  the  use  of  very  low 
frequency  transmissions^ ) ,   This  has  become  a  standard  pro- 
cedure throughout  the  world,  with  WWVL,  NBA  and  GBR  providing 
carriers  of  high  frequency  stability  for  scientific  use. 

The  procedure  for  obtaining  these  drift  rates  is  not 
complex,  but  because  of  the  small  errors  involved,  long  in- 
tervals of  measurement  are  necessary.  A  24  hour  measurement 
interval  is  recommended,  since  propagation  conditions  between 
receiver  and  transmitter  recur  at  24  hour  intervals.   In  ad- 
dition a  daylight  propagation  path  is  preferred  because  of 
the  more  stable  conditions.  Figures  11,  12,  13,  and  14  are 
typical  plots  of  drift  rates.  These  were  obtained  at  the 
U.  S.  Naval  Postgraduate  School.   A  tracking  receiver  with  a 
direct  reading  microseconds  dial  is  recommended  for  this  pur- 
pose, although  the  same  data  can  be  obtained  using  a  chart 
recorder.  At  24  hour  intervals  the  change  in  the  output  of 
the  phase  comparator  is  noted.   Since  this  output  is  cali- 
brated in  microseconds  it  can  be  used  directly  in  computa- 
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tionso   For  a  period  of  time  T,  the  error  iss 

error  ■    A  t 

(T)  (0„36) 

where  At  is  the  microseconds  change  during  T, 
When  several  day  have  been  plotted,  the  rate  of  drift  should 
begin  to  be  apparento   On  the  figures  of  drift  rates  referred 
to  above s  the  drift  is  shown  as  a  dotted  line,,   NPM  and  NAA 
did  not  appear  to  use  oscillators  of  sufficient  stability, 
and  drift  rates  could  not  be  determined  for  Lhese  stations „ 
Once  the  drift  rates  are  known,  corrections  can  be  predicted 
and  applied  to  the  navigational  calculations 0   To  find  the 
microseconds  correction 

correction  ■  (predicted  error)  (H)  (Q036)  (10-*-^)  where 
the  correction  is  in  microseconds  and  H  is  the  time  interval 
for  the  correction  (usually  the  elapsed  time  since  the  pre- 
vious f ix) o   As  an  example,  the  predicted  error  for  NBA  on 
21  January  1965  would  be  <=lo0  XIO"-1-^  and  a  correction  to  be 
applied  to  a  24  hour  time  period  would  be 

correction  -  (l.OxlCT10)  (24)  (Oo36)  (1010)  -  3„65  /I  sec 
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10c  Great  Circle  Calculations « 

The  microseconds  change  recorded  by  a  tracking  receiver, 
when  corrected  and  changed  to  distance s    represents  the  shift 
in  the  position  of  the  ship  along  the  great  circle  between 
the  ship  and  the  transmitter,.   This  change  is  applied  along 
the  great  circle  propagation  path  through  the  last  known 
position  of  the  ship0   If  a  circle  is  then  drawn  through  this 
corrected  pointy  with  the  center  of  the  circle  at  the  trans- 
mitter site,  the  circle  is  then  the  line  of  position  to  be 
used  for  the  new  fix0 

Charts  can  be  overprinted  with  circles  concentric  to 
the  VLF  stations  considered  useful  for  the  locality o   Figure 
15  is  a  portion  of  such  a  chart 0   In  this  case  the  circles 
are  33 °4  miles  (200^ sec)  apart,  with  the  distance  to  the 
station  noted  on  each  circle <>  When  a  new  ships  position  is 
desired,  the  lines  of  position  are  easily  found  by  drawing 
lines  parallel  to  the  circles  at  the  appropriate  distance 
from  the  last  known  position,, 

If  overprinted  charts  are  not  available,  the  great 
circle  bearing  and  distance  to  the  transmitter  can  be  found 
by  computation*.  Methods  for  computing  these  values  for  the 
spherical  earth  can  be  found  in  books  on  navigation  or  radio 
wave  propagat  ion ( 1 H  9 ) 0 

At  the  distance  where  the  VXF  signals  are  useful  the 
concentric  circles  approximate  straight  lines  for  small 
areaso   For  this  reason  the  line  of  position  can  be  drawn  as 
a  straight  line  perpendicular  to  the  computed  great  circle 
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FIG.      15 o      A   PORTION   OF    =  07". 

VLF  NAV]  ;      :on 


40 


patho  Unfortunately j  this  straight  line  approximation  will 
deviate  excessively  from  the  actual  concentric  circle  if  the 
new  position  is  more  than  about  50  miles  from  the  computed 
great  circle „   In  addition,  errors  may  result  from  the 
spherical  earth  calculation  since  the  earth  is  actually 
slightly  ellipticalo  With  this  form  of  relative  navigation 
such  errors  may  accumulate  over  a  period  of  days0 

In  order  to  keep  errors  at  a  minimum.,  many  of  the  great 
circle  calculations  are  necessaryc   These  can  become  tedious 
if  done  manually.,   Luckily  the  high  speed  digital  computer 
is  well  suited  to  these  calculations  *   These  calculations  can 
be  made  for  each  degree  of  latitude  and  longitude  and  the 
values  published  in  volumes  similar  to  other  tables  used  by 
navigators o   Errors  can  be  minimized  by  averaging  between 
the  given  values 0 
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lie   The  Ship  Navigation  Problem,, 

Many  changes  in  a  VLF  signal  that  are  apparent  and 
easily  corrected  at  a  fixed  receiver  location  are  not  so  ap- 
parent on  a  moving  ship  and  can  result  in  large  navigation 
errors o   Examples  of  these  areg 

lo  The  transmitter  leaves  the  air^  either  for  a  sched- 
uled maintenance  period  or  for  some  other  reason. 

2o  Cycle-slipping  occurs0 

3o  The  transmitter  master  oscillator  is  corrected. 

4o  The  transmitting  frequency  is  changed,, 

5o  Keying  is  changed  from  CW  to  FSKe 
Corrections  for  the  errors  caused  by  these  changes  are 
possible ,  and  methods  of  finding  these  corrections  are  de- 
scribed belowo 

The  VLF  stations  used  for  navigation  must  leave  the  air 
occasionally  for  routine  maintenance,,   During  these  periods 
the  tracking  receiver  on  the  ship  will  stop  tracking,  and 
when  the  station  resumes  transmitting  a  few  hours  later  the 
elapsed  microseconds  will  not  be  correct „   The  transmitter 
master  oscillator  is  not  stopped  during  these  maintenance 
periods  of  course ?  and  when  the  transmitter  comes  back  on 
the  air  the  phase  difference  between  the  local  oscillator 
and  the  received  signal  is  exactly  what  it  would  have  been 
if  the  transmitter  had  never  been  off •   At  a  fixed  receiver 
location  these  transmitter  outages  present  no  problem,. 
Figure  16  shows  an  outage  between  1330  and  1520  (simulated 
by  disconnecting  the  antenna )0   On  a  ship  the  problem  is  more 
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complex,  since  one  or  more  full  cycles  may  be  lost  during  the 
outage  because  of  the  movement  of  the  ship.   Consequently 
when  tracking  is  resumed  the  next  line  of  position  must  be 
changed  by  a  distance  equivalent  to  some  integral  number  of 
cycles.   If  the  length  of  the  outage  is  known,  the  number  of 
cycles  lost  can  be  computed  by  dividing  the  distance  the  ship 
has  traveled  away  or  toward  the  transmitter  by  the  length  of 
one  cycle.   If  the  length  of  the  outage  is  not  known,  or  if 
the  distance  traveled  cannot  be  determined,  the  corrected 
line  of  position  is  that  corrected  line  that  falls  nearest 
the  presumed  position  of  the  ship.   If  3  or  more  stations  are 
being  tracked  the  corrected  line  of  position  would  be  the 
line  that  falls  nearest  the  intersection  of  the  two  uncor- 
rected lines  (assuming  that  only  one  transmitter  left  the  air 
since  the  last  fix).  Table  VII  is  a  listing  of  some  wave- 
lengths in  microseconds,  kilometers  and  nautical  miles.   Val- 
ues for  frequencies  not  listed  in  Table  VII  can  be  found  by 
using  the  following  equations: 

t  =  1000 
f 

d  (miles)  =  t 


6.82 


d  (km)  =  t   =  d 

3.68    (miles) (1.852) 


where 

t  =  microseconds  per  wavelength 

f  =  frequency  in  kilocycles  per  second* 
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Frequency 

Wavelength 

f 
(kcs) 

t 
Usee) 

D 
(naut.   miles) 

D 
(km) 

16.0 

62.5 

9.16 

16.97 

17.3 

56.3 

3.25 

15.27 

13.6 

53«S 

7oS6 

14c55 

19»S 

50.5 

7o40 

13o70 

20.0 

50.0 

7o33 

13  c  57 

21.4 

46.75 

6. 34 

12.65 

24.0 

41o7 

6.11 

11.31 

26.1 

3S.3 

5o6l 

10,40 

TABLE  VII 


Wavelengths  of  Some  Commonly  Used  Frequencies 
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Frequency 

Wavelength 

f 
(kcs) 

* 

iusec) 

D 
(nauto  miles) 

D 
(km) 

I60O 

6205 

9ol6 

16,97 

17o3 

56o3 

So25 

15o27 

13o6 

53oS 

7o86 

14o55 

19oS 

50o5 

7o40 

13*70 

20oO 

50oO 

7o33 

13  c  57 

21.4 

46o75 

60S4 

12o65 

24o0 

41.7 

60II 

ilo31 

26ol 

3S.3 
1 

5o6l 

10o40 

TABLE  VII 


Wavelengths  of  Some  Commonly  Used  Frequencies 


45 


Cycle-slipping  is  another  possibility  for  erroru   This 
phenomenon  may  occur  occasionally  during  the  steeply  rising 
and  falling  portions  of  the  diurnal  phase  shift  trapezoido 
During  these  periods  the  diurnal  phase  shift  may  be  larger 
than  a  half  cycle  and  so  rapid  that  the  receiver  does  not 
track  the  proper  cycle  and  slips  a  full  cycle  to  the  nearer 
but  improper  tracking  point o   Another  cause  for  cycle-slip- 
ping may  be  phase  interference  between  signals  arriving  by 
both  great-circle  paths  or  from  multimodal  interference  dur- 
ing the  transition  from  night  to  day  or  day  to  night(24)0 
Figure  17  shows  cycle-slipping  on  the  track  of  NPM  (26ol  kcs) 
during  sunrise  at  the  UQ  S0  Naval  Postgraduate  School „   The 
dotted  line  shows  the  true  tracko   Correction  for  cycle- 
slipping  is  similar  to  the  correction  for  loss  of  signal 
since  an  integral  number  of  cycles,  normally  only  one?  have 
been  gained  or  losto   It  may  not  be  easy  to  spot  cycle- 
slipping  on  a  moving  ship,  so  a  correction  may  be  required 
if  the  line  of  position  seems  in  error  by  one  or  more  cycles <> 
Once  again 9   the  correction  will  be  more  obvious  if  many 
stations  are  being  trackedo 

Occasionally  it  may  be  necessary  to  readjust  the  trans- 
mitter master  oscillator  to  bring  it  within  tolerance 0  When 
this  is  dones  the  ships  record  will  be  in  error  by  the 
number  of  microseconds  equivalent  to  the  change 0   If  high 
quality  oscillators  are  used  at  the  transmitters^,  these 
changes  should  be  infrequent 0    In  addition  they  can  be 
scheduled  a  few  weeks  in  advance  and  ships  using  these  sta- 
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tions  for  navigation  may  be  notified  by  message „   It  might  be 
possible  to  correct  the  line  of  position  shortly  after  the 
change  is  made  by  adding  or  subtracting  the  microseconds 
equivalent 9  but  more  likely  the  correction  will  be  made  by 
forcing  the  line  of  position  through  the  first  fix  after  the 
oscillator  adjustment 0  The  oscillator  drift  rate  will  not  be 
affected  by  this  change 0 

The  Navy  operated  transmitters  may  change  frequency 0 
This  will  be  anticipated,  and  a  notification  of  this  change 
will  be  sento   Figure  l£  is  an  example  of  such  a  notice 0 
When  a  station  is  shifted  in  frequency  its  previous  tracking 
history  is  of  little  use  in  finding  a  future  line  of  posi- 
tion,,  A  new  fix  must  be  found  using  other  means ,  and  track- 
ing started  from  that  fix0   Since  the  stabilized  oscillator 
at  the  transmitter  is  not  altered  when  the  frequency  shift 
occurs,  the  drift  rate  at  the  new  frequency  should  remain 
constant o  When  many  stations  shift  frequency  as  in  Figure 
18,  the  new  fix  may  have  to  be  found  by  means  other  than 
VLF  transmissions,, 

When  the  transmitted  signal  is  shifted  from  CW  to  FSK 
type  keying^  it  may  or  may  not  be  usable  for  navigation,. 
Measurements  at  the  Uc  So  Naval  Postgraduate  School  have  in- 
dicated that  NAA  can  not  be  tracked  while  transmitting  FSK9 
while  NPM  can  be  used  for  navigation  even  when  transmitting 
FSK,  provided  a  10  microsecond  correction  is  added „  The  rea= 
son  for  this  may  be  that  NAA  uses  the  wider  frequency  shift, 
or  NPM  may  have  one  side  of  the  frequency  shift  stabilized 
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U.  S.  NAVAL  OBSERVATORY 
WASHINGTON,  D.   C.     20390 


TIME  SERVICE  ANNOUNCEMENT  21  May  1964 

Change 8  in  the  VLF  Transmissions  of 


NBA,  NLK/NPG.   and  NAA 


1.  The  frequency  of  NBA,  Balboa,   Canal  Zone,  will  be 
changed  from  18.0  kc/s  to  Z4.0  kc/s  on  8  June  1964. 

2.  The  frequency  of  NLK/NPG,  Jim  Creek,  Washington,  will 
be  changed  from  24.0  kc/s  to  18.6  kc/s  on  1  June  1964. 

3.  The  frequency  of  NAA,  Cutler,  Maine,  will  be  changed  to 
17.8  kc/s  with  FSK  when  transmission  is  resumed  on  5  June  1964, 
and  will  not  be  useable  for  frequency  calibration* 

4.  The  frequency  of  NSS,  Annapolis,  Maryland,  continues  on 
21.4  kc/s. 


T.  S.  BASKETT 
Superintendent 


FIGo      16\      A  FREQUENCY  CHANGE  NOTICE 
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to  the  precision  frequency  standard,.   If  the  use  of  VLF 
transmissions  for  navigation  becomes  common,  all  FSK  trans- 
missions will  probably  be  made  usable o 

On  a  ship,  the  major  problems  in  tracking  VLF  stations 
for  navigation  occur  because  the  ship  is  moving ,  and  most  of 
the  corrections  necessary  require  a  good  fix0   For  this  rea- 
son it  will  be  useful  to  track  as  many  VLF  stations  as  pos- 
sible in  order  to  have  a  good  fix  available  when  there  is 
some  change  in  a  transmission  during  the  run  between  fixes „ 
Another  problem  is  the  high  frequencies  presently  used  by 
most  Navy  transmitters..   The  lower  frequencies  are  the  most 
useful  for  navigation  because  there  is  less  probability  of 
an  error  on  one  cycle  when  a  wavelength  is  long  with  respect 
to  the  possible  accumulated  errors „   As  indicated  in  Table 
VII,  one  cycle  for  NPM  at  26»1  kcs.  occurs  every  5°6l  miles, 
and  this  frequency  could  very  likely  be  too  high  to  be  usedo 

An  important  feature  of  VLF  navigation  that  must  be 
recognized  is  that  an  error  in  a  position  will  not  correct 
itself  at  a  subsequent  fixo   Once  a  ships  position  is  found, 
the  lines  of  position  are  carried  forward  from  that  point, 
and  if  this  fix  were  in  error  this  error  would  be  carried  on 
to  the  next  fix  with  no  change »   Because  of  this,  it  will  be 
necessary  to  periodically  check  the  snipes  position  as  found 
from  VLF  by  other  independent  means* 

One  of  the  advantages  of  the  VLF  navigation  system  is 
the  possibility  of  an  automated^,  accurate ,  all  weather  navi- 
gation system  using  this  scheme «   It  should  be  possible  to 
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continuously  plot  the  ships  position  on  standard  navigation 
charts.   Although  the  instrumentation  of  such  a  system  was 
not  investigated,  it  does  appear  possible  to  construct  such 
a  device,,   Since  corrections  for  most  common  sources  of  error 
are  similar  (searching  through  integral  wavelengths  for  the 
proper  line  of  position)  a  program  to  compensate  for  these 
can  be  developed .  An  operator  would  probably  have  to  insert 
corrections  for  other  possible  errors ,  and  the  output  of  the 
system  would  have  to  be  checked  against  a  known  position 
periodically  to  eliminate  accumulated  error., 

A  particularly  useful  application  of  the  VLF  navigation 
system  is  its  use  on  nearly  stationary  ships,  such  as  a  ship 
on  a  semi-permanent  ocean  station,,  When  used  at  a  fixed 
location  any  errors  contained  in  a  received  signal  can  easily 
be  spotted  and  corrected,,   In  addition,  24  hour  navigation  is 
possible  since  the  previous  days  track. of  a  VLF  signal  will 
provide  accurate  corrections  for  the  diurnal  phase  shift 
during  the  following  24  hour  period 0 
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12.  A  Sample  Problem,, 

Finding  a  ships  position  from  the  VLF  transmission  is 
not  difficult,  provided  overprinted  charts  or  tables  of  VLF 
bearings  and  distances  are  available.  Figure  19  is  an  ex- 
ample of  a  page  from  such  a  table. 

A  typical  problem  is  shown  on  a  portion  of  a  position 
plotting  sheet  in  Figure  20.  The  ships  position  is  known  at 
0300.  This  position  could  have  been  obtained  from  a  star 
fix,  Loran,  or  a  previous  VLF  fix.  The  phase  difference 
reading  for  the  VLF  stations  of  interest  are  also  known. 
This  information  is  obtained  either  by  recording  the  micro- 
seconds dial  reading  on  the  VLF  receiver  at  0800,  or  by 
finding  the  readings  at  this  time  on  the  charts  recording  the 
output  of  the  phase  comparators.   It  is  desired  to  fix  the 
position  of  the  ship  at  1200.  Table  VIII  shows  the  bearing 
and  ranges  for  the  VLF  stations.  These  were  obtained  from 
tables  similar  to  Figure  20.  The  values  are  for  the  in- 
teger latitude  and  longitude  nearest  the  known  ships  position 
at  0300,  in  this  case  123. 0  W,  33.0  N,  From  Table  VIII  we 
that  NLK,  NPM,  and  WWVL  are  within  the  near  zone  of  4000 
kilometers  (2160  N  miles),  and  signals  from  these  stations 
will  not  be  used. 

It  will  be  necessary  to  compute  the  expected  diurnal 
phase  shift  for  the  remaining  stations.  The  times  of  sunrise 
and  sunset  at  each  of  the  usable  stations  must  be  known,  and 
transferred  to  local  time.  These  values  may  be  available 
from  previous  computed  tables,  or  they  can  be  computed  on 
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TRANSMITTER  STATION-GBR 


RECFIVER 

POSITION 

BEARING 

DISTANCE 

LAT 

LONG 

(NMILES) 

(KM) 

10 .  on 

129.  no 

29 

.72 

6208.58 

11498.30 

11.  00 

129.  00 

29 

.60 

6156.73 

11402.27 

12,00 

129. 00 

29 

,5n 

61 04.82 

11306.13 

13.00 

129.  00 

29 

.40 

6052.85 

11209.88 

14.  00 

129 . no 

29 

.31 

6000.83 

11113.54 

15.00 

129.00 

29 

.23 

5948.76 

11017.11 

16.00 

129.  00 

29 

.16 

5896.65 

10920.60 

17.00 

129. 00 

29 

.09 

5844.50 

10824.01 

18.  00 

129. 00 

29 

.03 

5792.31 

10727.36 

19.00 

129.00 

28 

.98 

5740  .  09 

10630.64 

20.  00 

129. 00 

28 

.94 

5687.83 

10533.87 

21.  00 

129. 00 

28 

.90 

5635.56 

10437.05 

22.  on 

129.  no 

28 

.88 

5583.26 

10340.19 

23.  00 

129.no 

28 

.86 

5530.94 

10243.29 

24.00 

129.00 

28 

.84 

5478.60 

10146.37 

25.00 

129.  00 

28 

.84 

5426.25 

10049.42 

26.00 

129.  00 

28 

.84 

5373.90 

9952.46 

27.00 

129.00 

28 

.85 

5321.54 

9855.48' 

28.00 

129.  00 

28 

.86 

5269.17 

9758.51 

29.00 

129. 00 

28 

.89 

5216.81 

9661.54 

30.  00 

129.  00 

28 

.92 

5164.46 

9564.57 

31.00 

129.  no 

28 

.96 

5112.11 

9467.63 

32.  00 

129.no 

29 

.no 

5059.78 

9370.71 

33.  00 

129.  00 

29 

06 

5007.46 

9273.82 

34.  on 

129.  no 

29 

.12 

4955.16 

9176.96 

35.00 

129.  00 

29 

19 

4902.89 

9080.16 

36.  0  0 

129.  00 

29 

26 

4850.65 

8983.40 

37.  00 

129. no 

29 

.35 

4798.44 

8886.70 

38.  0  0 

129.  no 

29 

44 

4746.26 

8790.07 

39.  00 

129.  00 

29 

54 

•  4694.12 

8693.51 

40.  on 

129. no 

29 

65 

4642.03 

8597.  04 

41.00 

129. no 

29. 

77 

4589,99 

8500.66 

42.  00 

129. no 

29. 

89 

45  38. 0  0 

8404.37 

43.  on 

129.00 

30. 

02 

4486.  07 

8308.19 

44.  00 

129. no 

30 

17 

4434.20 

8212.13 

45.  00 

129. 00 

30. 

32 

4382.4Q 

8116.20 

46.  00 

129.  00 

30. 

4fl 

4330.67 

8020.40 

47.  00 

129.00 

30. 

65 

4279.02 

7924.74 

48.  00 

129.00 

30. 

83 

4227.45 

7829.24 

49.  00 

129.00 

31. 

02 

4175.97 

7733.90 

50.  on 

129.  no 

31. 

22 

4124.59 

7638.74 

51.  00 

129. 00 

31. 

43 

4  073.31 

7543,7.7 

52.  00 

129.  00 

31. 

64 

4022.14 

7449.00 

53.  00 

129.  no 

31. 

87 

3971.08 

7354.44 

54.  on 

129.  00 

32. 

12 

3920.14 

7260.10 

55.  00 

129.  00 

32. 

37 

3869.33 

7166.  00 

56.  on 

129.00 

32. 

63 

3818.66 

7072.16 

57.  on 

129. no 

32. 

91 

3768.13 

6978.58 

58.  00 

129. no 

33. 

20 

3717.75 

6885.28 

59.00 

129.00 

33. 

50 

3667,54 

6792.28 

60.  00 

129. 00 

33. 

81 

3617.49 

6699.59 

FIG.   19.   A  PAGE  FROM  A  TABLE  OF  GREAT  CIRCLE  RANGES 
AND  BEARINGS 
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board  ship0   The  times  of  sunrise  and  sunset  locally  must  be 
calculatedo   The  times  of  sunrise  and  sunset  calculated  and 
listed  in  the  tables  are  not  those  at  sea  level,  but  those  at 
an  altitude  of  70  kilometers  where  the  lower  reflecting  layer 
of  the  ionosphere  forms o   Values  of  kv  for  the  local  area 
must  also  be  known  so  the  A  t  may  be  calculatedo 


Station 

Bearing 
(degrees) 

Distance 
(N  miles) 

NLK 

21o4# 

666069 

NPM 

246o 04 

1347oU 

NSS 

71o9S 

23 94 0 26 

NAA 

6lo31 

2703 0 03 

NBA 

110o70 

3129.30 

GBR 

30o02 

4670 0 90 

WWVL 

74ol9 

IO760 47 

TABLE  VIII 


Bearing  and  Distance  of  VLF  Stations 
from  128o0W8  3S0O  N 


Table  IX  shows  the  times  of  sunrise  and  sunset  for  the 
desired  VLF  stations  for  a  typical  day0 
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VLF 
Station 

GMT 
Sunrise      Sunset 

NBA 

1052 

2351 

NSS 

1035 

2343 

NAA 

0956 

2307 

GBR 

0523 

1346 

TABLE  IX 


Times  of  Sunrise  and  Sunset  for  a 
Typical  Day 


Local  sunrise  and  sunset  are  computed,  corrected  for  altitude 
and  found  to  be  052S  and  1&36  respectively,  local  time..  The 
tabulated  values  of  sunrise  and  sunset  from  Table  IX  are  cor- 
rected to  local  time  by  adding  the  time  equivalent  of  the 
differences  in  longitude. 

123 oOW  Longitude  =  8hr.  32  min. 

local  times  of  sunrise  and  sunset  for 

the  VLF  stations 


sunrise 

sunset 

NBA 

0220 

1519 

NSS 

0203 

1511 

NAA 

0124 

1435 

GBR 

2056 
56 

1014 

Values  of  ky  are  found  from  Table  II  for  NSS,  NBA,  and  GBR 
The  value  of  £  t  is  found  from  this  and  the  values  of  dis- 
tances in  Table  VII,  £t  =  (ky  D)(loS6)0 


VLF  Station 

Freq 
(kcs) 

(xlS3) 

At 

U^sec 

NBA 

24o0 

So20 

47.7 

NSS 

21o4 

7.65 

34.0 

NAA 

17oS 

7o22 

36c4 

GBR 

I6e0 

7.2S 

63 .3 

The  values  of  ky  used  are  the  best  available  and  will  serve 
as  examples o 

From  the  values  found  for  sunrise,  sunset  and  £  t,  the 
trapezoidal  pattern  expected  for  the  diurnal  shift  can  be 
constructed o   These  are  shown  in  Figure  21 „ 

From  Figure  21  it  is  clearly  seen  that  during  the  period 
0300-1200  local  time,  the  propagation  paths  for  NBA,  NSS  and 
NAA  are  completely  in  daylight  and  no  correction  for  diurnal 
shift  is  necessary  for  these  stations »  The  GBR  propagation 
path  is  in  daylight  at  0800  but  is  partially  in  the  dark  at 
1200c   The  correction  to  the  microseconds  reading  at  1200  can 
be  found  from  Figure  21,  and  is  13  microseconds 0 

The  final  correction  to  the  microseconds  difference 
reading  is  that  correction  for  drift  rate0   These  rates  are 
obtained  as  explained  in  section  9  from  graphs  similar  to 
Figures  13 ,  14>  15  and  16 •   Some  typical  drift  rates  might  be 
as  follows,  where  drift  rate  is  in  microseconds  per  24  hours e 
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FIG.   21.   CONSTRUCTED  TRAPEZOIDAL  DIURNAL  SHIFT  PATTERNS 
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Station     Drift  Rate     Drift  in  4  Hour 
C#sec/24hr)        (^sec) 
NBA  »10  -1.67 

NSS  -8  -1.33 

NAA  -10  -1.67 

GBR  -9  -1.50 

The  corrections  are  now  ready  to  be  applied  to  the 
elapsed  microseconds  and  the  various  lines  of  position  con- 
structed. The  elapsed  microseconds  between  0800  and  1200 
might  be  found  to  be 

VLF  Station        Elapsed  Microseconds 
NBA  +334.7 

NSS  +194.1 

NAA  +143.8 

GBR  -  64.6 

The  corrections  are  applied  and  the  change  in  ships  position 
along  the  great  circle  propagation  path  computed. 

NBA        NSS       NAA       GBR 

Elapsed  sec           +334.7     +194.1    +143. 8     -64.6 
Diurnal  correction     0     0        0     +13 .0 


Drift  rate 
Corrected  sec 


+334.7     +194.1    +143.8     -51.6 
-1.7      -1.3      -1.7     -1.5 
+333.0     +192.8    +142.1     -53.1 


N  Miles  Change        +53.8     +  31.2    +23.0     -  9.4 
These  nautical  mile  changes  are  applied  along  the  great 
circle  path,  away  from  the  transmitter  for  a  positive  change 
and  toward  the  transmitter  for  a  negative  change.   It  is  use- 
ful,  but  not  essential,  to  have  a  rough  estimate  of  the  new 
ships  position  from  a  dead-reckoning  plot  so  that  changes  in 
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the  great  circle  bearings  can  be  averaged  to  reduce  error. 
In  this  sample  problem  the  new  ships  position  is  known  to  be 
near  129«0W,  39oON,  and  the  great  circle  ranges  and  bearings 
are  obtained  from  the  tabulated  values.   These  are  listed  in 
Table  X0  Comparing  these  values  with  those  in  Table  VIII  it 
is  seen  that  there  is  little  difference  in  bearing  for  those 
stations  beyond  4000  kilometers.   The  difference  in  bearing 
caused  by  averaging  will  certainly  be  less  than  the  con- 
struction error  on  Figure  20.  These  averages  are  found  how- 
ever and  ares 

Station      Average  Great  Circle  Bearing 

NBA  110.64 

NSS  72.32 

NAA  61.62 

GBR  29. 73 

These  bearings  are  plotted  and  the  computed  changes  in  miles 
are  measured  along  these  lines  in  the  proper  direction.   The 
lines  of  position  are  then  drawn  perpendicular  to  the  great 
circle  lines  at  this  point.   The  intersection  of  these  lines 
of  position  results  in  a  fix,  as  shown  in  Figure  20. 
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Station 

Bearing 
(degrees) 

.  . ,.    .  .     _ 

Distance 
(N  miles) 

NLK 

26o  74 

631o30 

NPM 

243 o32 

1830.30 

NSS 

72o77 

2420o99 

NAA 

6lo93 

2721o01 

NBA 

110o5B 

3194oS6 

GBR 

29o54 

4694c 12 

WWVL 

77  =  12 

1107 02S 

TABLE  X 


Bearing  and  Distance  of  VLF  Stations 
from  129oOW„  39«ON 


The  procedure  for  finding  the  ships  position  in  an  area 
where  overprinted  charts  are  available  is  similar  except  that 
tables  of  bearing  and  distance  to  the  VXF  stations  are  not 
necessary  <>   This  information  is  contained  on  the  chart  in  the 
form  of  lines  of  equal  distance  from  the  stations 0   The  nau<= 
tical  miles  change  is  computed  as  before ,  and  this  change  is 
measured  along  a  line  through  the  last  ships  position  perpen- 
dicular to  the  lines  of  equal  distance  for  the  VXF  station,, 
The  line  of  position  is  then  a  line  parallel  to  these  lines 
of  equal  distance • 
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13.   Comparison  With  Other  Electronic  Navigation  Systems. 

It  appears  pertinent  to  examine  how  the  VLF  navigation 
system  of  this  paper  compares  generally  with  other  electronic 
navigation  systems.   Some  other  systems  presently  in  use  are 
LORAN  A,  LORAN  C,  DECCA  and  Radux-Omega. 

LORAN  A  operates  at  frequencies  near  2  megacycles  and  is 
useful  at  ranges  out  to  about  650  nautical  miles  during  day- 
light and  about  twice  that  at  night.   It  is  a  pulsed  system 
utilizing  a  master  and  a  slave  station.   A  special  receiver 
measures  the  difference  between  the  time  of  arrival  at  the 
ship  of  the  pulses  from  the  master  and  slave  signals.   Spe- 
cial charts  are  published,  giving  lines  of  constant  time 
difference  for  each  pair  of  stations  within  range.   The  posi- 
tion of  the  ship  is  found  at  the  intersection  of  the  time 
difference  lines.   Many  pairs  of  stations  are  required  to 
cover  an  ocean  area  because  of  the  relatively  short  useful 
range  of  the  transmitted  signal.   Propagation  inconsistencies 
can  cause  large  errors  at  distances  where  reception  of  the 
pulses  is  marginal. 

LORAN  C  is  also  a  pulsed  system,  but  operates  at  about 
100  kcs„,  and  consequently  has  longer  range  with  more  stable 
propagation  conditions.   In  addition  more  accurate  comparison 
of  the  two  received  pulses  is  possible  since  the  LORAN  C  re- 
ceiver compares  the  time  between  separate  and  corresponding 
oscillations  of  the  master  and  slave  stations.   This  is 
called  cycle-matching.   As  with  LORAN  A,  a  fix  is  found  at 
the  intersection  of  time  difference  lines  on  special  LORAN  C 
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charts.   LORAN  C  is  considered  useful  up  to  a  2000  mile  range 
with  good  accuracy0 

The  DECCA  system  also  operates  at  low  frequency,  but 
does  not  use  pulses c   In  this  system  a  master  station  is  sur- 
rounded by  three  slave  stations  at  distances  of  60  to  120 
nautical  miles »  These  stations  operate  at  frequencies  that 
are  all  different,  but  are  all  harmonics  of  some  lower  funda- 
mental frequencyc   The  transmitter  frequencies  may  lie  in  the 
range  of  70  to  120  kcs<,   The  four  transmitters  broadcast  a 
continuous  wave,  with  all  slave  stations  related  in  phase  to 
the  master  station0   Special  equipment  aboard  the  ship  com- 
pares the  phase  between  the  four  stations  and  gives  an  output 
on  special  recorders,  called  decometerso   A  fix  is  obtained 
by  transferring  the  decometer  readings  to  special  charts  with 
lines  corresponding  to  various  readings „   A  ships  position  is 
found  at  the  intersection  of  the  indicated  lines*   In  spite 
of  the  low  operating  frequency j  this  is  considered  to  be  a 
short  range  system,  with  a  maximum  range  of  240  nautical 
miles o   Accuracy  varies  from  less  than  one  quarter  mile  dur- 
ing full  daylight  at  close  range  to  &   miles  during  dusk  at 
maximum  range ( H ) • 

The  Radux-Omega  navigation  system  operates  at  VLF  fre- 
quencies, and  is  a  long  range  system  using  phase  comparison 
between  signals  transmitted  by  a  master  and  a  slave  station,, 
The  slave  station  is  phase  locked  to  the  signal  transmitted 
by  the  master  station*   The  master  and  slave  stations  trans- 
mit alternately  on  the  same  frequency,  and  equipment  on  the 
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ship  measures  the  phase  difference  between  the  two  received 
signals*   Phase  ambiguities  occur  in  lanes  separated  by  one- 
half  the  wavelength  of  the  operating  frequency*   Lane  iden- 
tification is  accomplished  by  alternate  transmission  by  both 
master  and  slave  on  two  frequencies 0   A  ships  position  is  ob- 
tained at  the  crossing  point  of  lines  of  equal  phase  differ- 
ence.  Charts  are  printed  containing  the  lines  of  equal  phase 
difference  for  each  pair  of  stations0   Ranges  in  excess  of 
5000  nautical  miles  are  expected,  and  a  network  of  eight  sta- 
tions can  provide  15  sets  of  lines  of  position  by  phase  re- 
lations between  the  various  stations^^)  o 

The  method  of  navigation  discussed  in  this  paper  differs 
from  the  systems  previously  discussed „   This  system  can  be 
used,  for  example,  without  specially  prepared  charts  by  uti- 
lizing  tables  of  previously  computed  great  circle  calcula- 
tions 0      In  addition,  the  shore  stations  utilized  are  pres- 
ently in  operation  for  purposes  of  communications «   As  with 
Radux-Omega  the  received  signal  must  be  corrected  for  changes 
along  the  propagation  path,  but  in  this  system  the  correction 
can  be  found  from  the  previous  days  recordo   A  disadvantage 
is  the  lane  identification  problem,  but  this  can  be  simpli- 
fied by  utilizing  lower  operating  frequencies  to  widen  the 
lanes,  and  by  using  as  many  stations  as  possible  for  ob= 
taining  the  fix<>   The  lane  identification  problem  can  be 
further  simplified  by  using  other  methods  of  navigation  to 
find  the  ships  position  within  a  half  lane  widtho   The  ship- 
board equipment  necessary  for  this  system  is  also  simple  when 
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compared  to  the  equipment  necessary  for  most  existing  sys- 
tems. 

The  accuracy  of  this  system  should  be  close  to  that  pre- 
dicted for  the  Omega  systerrp-^).   In  the  Omega  system,  a  24 
hour  RMS  accuracy  of  0.5  mile  during  the  day  and  1.2  miles  at 
night  is  predicted o   For  24  hour  use  of  the  Omega  system, 
tables  of  corrections  for  the  diurnal  phase  shift  must  be 
available. 

The  method  of  navigation  discussed  in  this  paper  may  not 
be  usable  by  Naval  ships  that  change  course  often  and  spend 
long  periods  at  sea.   In  its  present  form  is  should  be  par- 
ticularly useful  for  ships  that  desire  to  remain  at  a  fixed 
ocean  station  or  by  oceanographic  vessels. 
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compared  to  the  equipment  necessary  for  most  existing  sys- 
teraso 

The  accuracy  of  this  system  should  be  close  to  that  pre- 
dicted for  the  Omega  system'12) o   In  the  Omega  system,  a  2k 
hour  RMS  accuracy  of  0„5  mile  during  the  day  and  1<,2  miles  at 
night  is  predicted o   For  24  hour  use  of  the  Omega  system^ 
tables  of  corrections  for  the  diurnal  phase  shift  must  be 
available,, 

The  method  of  navigation  discussed  in  this  paper  may  not 
be  usable  by  Naval  ships  that  change  course  often  and  spend 
long  periods  at  sea0   In  its  present  form  is  should  be  par- 
ticularly useful  for  ships  that  desire  to  remain  at  a  fixed 
ocean  station  or  by  oceanographic  vessels,, 
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14.   Comments  and  Recommendations, 

The  primary  test  of  a  navigation  system  is  how  well  a 
ships  navigator  can  use  it  for  navigation   Most  of  the  data 
used  in  this  paper  were  obtained  at  the  U0  S0  Naval  Postgrad- 
uate School  at  Monterey,  California c   It  is  not  difficult  to 
navigate  a  building,  which  is  one  way  to  say  the  author  was 
not  able  to  test  the  system  properly0   As  a  consequence,  be- 
fore any  consideration  can  be  given  to  the  value  of  this  sys- 
tem, considerable  testing  under  shipboard  conditions  will  be 
necessary.   For  instance,  how  does  the  roll  and  pitch  of  a 
ship  using  a  vertical  antenna  for  reception  effect  the  recep- 
tion of  a  VLF  signal  of  very  low  signal  strength,,   This 
testing  seems  appropriate  and  is  recommended.. 

The  data  on  drift  rates  v/hich  is  plotted  as  Figures  11, 
12,  13  and  14  was  obtained  using  a  Rubidium  Vapor  Frequency 
standard,  as  noted  in  Chapter  7o   The  spread  in  data  points, 
although  generally  within  the  frequency  tolerance  of  the  VLF 
stations,  was  greater  than  expected <,   It  may  be  that  some  of 
the  spread  was  caused  by  the  local  standard  and  because  of 
this  a  good  crystal  oscillator  would  be  preferred  for  navi- 
gation purposes .   A  crystal  oscillator  of  sufficient  stabil- 
ity  was  not  available  during  periods  data  was  taken ,   When 
such  an  oscillator  is  available  similar  data  should  be  taken 
using  the  crystal  as  the  primary  standard „   If  possible  both 
the  crystal  oscillator  and  the  rubidium  standard  should  be 
used  to  obtain  drift  rates  simultaneously,,   Two  tracking  re- 
ceivers will  be  necessary  for  this  procedure 0 
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The  author  had  hoped  to  compare  antennas  used  for  the 
VLF  tracking  receiver.   The  antenna  used  was  a  Textran  loop, 
oriented  to  receive  signals  from  the  North-East.   A  loop  is 
not  satisfactory  for  use  on  a  ship  and  a  comparison  between 
the  loop  and  a  good  quality  whip  would  be  useful  when  a  whip 
becomes  available.  Unfortunately  the  noise  level  at  the  an- 
tenna location  used  at  the  U.  S.  Naval  Postgraduate  School  is 
high,  and  a  whip  may  not  be  usable  there. 

It  was  concluded  that  signals  within  the  multi-mode  zone 
were  not  usable  because  of  the  interaction  between  modes. 
This  is  unfortunate  since  the  VLF  signals  are  of  high  signal 
strength  in  this  zone.   In  particular,  if  VLF  signals  were 
usable  for  navigation  in  the  multi-mode  zone,  transmitters 
with  low  radiated  power  could  be  established  for  navigation 
use.   Since  efficient  antenna  systems  are  the  major  expense 
of  a  VLF  transmitting  station,  many  low  power  transmitters 
with  inefficient  antennas  can  be  established  at  relatively 
low  cost.  Reception  at  these  low  level  signals  is  not  diffi- 
cult with  the  narrow  bandwidth  tracking  receiver.  A  trans- 
mitter radiating  10  watts  at  a  frequency  of  15  kilocycles 
should  be  useful  for  navigation  at  ranges  in  excess  of  2000 
miles,  and  1  watt  of  radiated  power  may  be  used  at  distances 
beyond  1000  miles.  Even  0.1  watt  might  be  useful  at  consid- 
erable distances.   Of  course  it  may  be  difficult  to  receive 
signals  of  such  low  levels  unless  they  are  known  to  be 
present  at  the  exact  frequency.  The  multi-mode  zone  should 
be  investigated  to  determine  how  it  might  be  used  for  navi- 
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gation,  possibly  with  such  low  level  signals. 

Another  aspect  of  this  VLF  navigation  system  that  is 
deserving  of  further  investigation  is  the  development  of  an 
automated  navigation  system  using  VLF  phase  comparison.   This 
method  seems  particularly  adaptable  to  automatic  plotting 
since  phase  variations  transfer  linearly  to  distance  changes. 
Other  VXF  and  LF  systems,  as  noted  in  Chapter  13,  use  phase 
comparisons  between  master  and  slave  stations,  and  conse- 
quently the  lines  of  equal  phase  difference  are  hyperbolas. 

In  conclusion,  the  author  wishes  to  thank  the  many  in- 
dividuals who  gave  valuable  advice  and  counsel  on  subjects 
relating  to  this  paper.   Particular  appreciation  is  given  to 
Mr.  J.  H.  Stanbrough  Jr.  of  the  Woods  Hole  Oceanographic 
Institution  who  made  available  many  of  the  charts  and  tables 
used  in  the  1963  Indian  Ocean  cruise  of  the  ATLANTIS  II. 
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